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Abstract: We describe the lowest Landau level of a quantum electron star in AdS4. 
In the presence of a suitably strong magnetic field, the dynamics of fermions in the bulk 
is effectively reduced from four to two dimensions. These two-dimensional fermions can 
subsequently be treated using the techniques of bosonization and the difficult many- 
body problem of building a gravitating, charged quantum star is reduced to solving the 
sine-Gordon model coupled to a gauge field and a metric. The kinks of the sine-Gordon 
model provide the holographic dual of the lowest Landau levels of the strongly-coupled 
d — 2 + 1-dimensional boundary field theory. The system exhibits order one oscillations 
in the magnetic susceptibility, now arising as a classical effect in the bulk. Moreover, 
as the chemical potential is varied, we find jumps in the charge density, oscillations in 
the fractionalised charge density and plateaux in the cohesive charge density. 
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1. Introduction 

If you want to construct a gravitating star out of quantum fields, it's a lot easier to 
work with bosons than fermions. For bosons, you can place a macroscopic number of 
particles in the same state and, correspondingly, you need only find a single solution 
of the coupled Einstein-scalar equations. For fermions, the exclusion principle means 
this isn't possible: instead you need to solve the coupled Einstein-Dirac equations 10 23 
times. 



1 



Of course, to describe stars in the night sky, a fortuitous separation of scales means 
that you don't need to start with quantum fields. The Compton wavelength of the 
particles is much (much) less than the curvature of space. This ensures that you can 
treat the fermions as a fluid of particles, using standard statistical mechanics tech- 
niques to compute their equation of state and stress tensor which subsequently feed 
into the Einstein equations. However, in situations where the Compton wavelength of 
the particles is comparable to the curvature of space, you are left with a formidable 
many-body problem: the construction of a genuinely quantum star. 

The problem of constructing quantum stars has recently arisen in the search for 
holographic models describing a finite density of strongly interacting fermions. Initial 
attempts to understand fermions at finite density focussed on probe Dirac fields in an 
AdS black hole background [1, 2, 3, 4]. This class of models was shown to encompass 
a large variety of Fermi liquids, non-Fermi liquids and marginal Fermi liquids. See [5] 
for a detailed review. 

While the results of the probe calculations are extremely encouraging, they leave a 
number of unanswered questions. Most notably, the charge density in the boundary 
theory is not carried by the fermions themselves, but is rather contained in 'fraction- 
alised' charge, hidden behind the horizon of the black hole. Moreover, the horizon 
suffers from a non-zero entropy density at zero temperature, in violation of the third 
law of thermodynamics and strongly suggesting that the black hole is not the true 
ground state of the system. Indeed, one important instability was identified in [6]: 
the fermions want to develop a charge density outside the horizon. This charge den- 
sity is sufficient to cause appreciable backreaction, both screening the electric field of 
the black hole and, ultimately, replacing the horizon with a different geometry with 
vanishing entropy. The resulting configuration is referred to as an electron star [7]. 

The phenomenology of electron stars was subsequently explored in a number of pa- 
pers [8, 9, 10, 11]. However, to avoid the complexities describe above, all these papers 
work in a regime in which the mass, m, of the fermion is large compared to L, the 
AdS radius: mL 3> 1. This provides the necessary separation of scales to model the 
equation of the state of the star using a simple fluid description. While this makes the 
calculations computationally tractable, it is not without consequence. For example, the 
electron star exhibits a large number of densely packed Fermi surfaces, corresponding 
to the radial decomposition of the fermion field in the bulk [9, 12]. The existence of 
many Fermi surfaces is not particularly desirable and can be traced directly to the need 
to work in mL ^> 1 approximation. 



2 



More recently, progress has been made in constructing a truly quantum electron star 
in a regime where mL ~ 1 [13, 14]. If you're only interested in the geometry rather than 
the fermionic wavef unctions, you don't really need to solve the equations 10 23 times; 
you only need to calculate the determinant of the Dirac operator in a state which is 
occupied by 10 23 fermions. However, this is still a technically challenging subject and 
is not yet as well developed as the earlier approach. 

The purpose of this paper is to point out that there is a regime in which the full 
quantum electron star in AdS4 can be constructed with relative ease. This is the regime 
of large magnetic field, B. By this we mean B > fi 2 , where /i is the chemical potential 
of the boundary field theory. 

There is a simple, intuitive explanation for why electron stars are particularly simple 
in this limit. In the presence of a large magnetic field, the bulk electrons are forced 
into their lowest Landau level. This results in a kind of dimensional reduction of their 
dynamics: the electrons are trapped in the spatial directions parallel to the boundary 
and only free to move in the radial direction of AdS. In other words, the fermions 
become effectively two dimensional. 

But fermions in two dimensions are special. They are equivalent to bosons. This 
means that the problem involving a finite density of quantum fermions can be mapped 
to a classical problem of a scalar field, coupled to a gauge field and metric. And, as we 
mentioned in the opening paragraph, this is much simpler to solve. Indeed, at heart 
the problem of constructing a quantum electron star becomes morally equivalent to the 
problem of screening an electric field by fermions in two dimensions [15]: we will find 
that we need to solve a variant of the sine-Gordon model in AdS-like geometries. This 
approach also provides us with the analog of the fermion wavefunctions, providing an 
intuitive visualisation of the distribution of fermionic charge in the bulk. 

The boundary dual of our holographic model is a 2 + 1 -dimensional, strongly-coupled, 
field theory. The theory is studied at finite density, with a chemical potential /x, and 
in the presence of a magnetic field B. We are restricting ourselves to the limit 

B > /i 2 

If the boundary theory consisted of free fermions, this requirement would force us 
to the lowest Landau level. We will find that something similar holds for our strongly 
interacting theory. We are now dealing with large N gauge theory and, correspondingly, 
the gauge invariant fermionic operators have a large number of "lowest Landau levels" , 
corresponding to the radial harmonics of the fermionic field in the bulk. We will refer 
to these different Landau levels as different "bands" . 
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Our construction of electron stars in the lowest Landau level allows us to exhibit 
a number of quantum phenomena using purely classical bulk physics. These include 
discrete jumps in the charge density as the chemical potential is varied and de Haas 
van Alphen oscillations as the magnetic field is varied. 

The organisation of the paper is as follows. Section 2 is devoted to the bulk. We study 
the decomposition of fermions in a magnetic field and explain how one can bosonize 
those fermions in the lowest Landau level. We further review some well known aspects 
of bosonization, including the existence of kinks and their relationship to the anomaly, 
both of which gain a slight twist in the present context. In Section 3, we solve our 
bosonized equations to construct quantum electron stars in both hard wall and black 
hole geometries. We subsequently use these solutions to explore various phenomena 
in the d — 2 + 1 dimensional boundary theory, such as oscillations in the magnetic 
susceptibility and the interplay between filling fractions and the fractionalisation of 
charge due to the black hole horizon. A number of technical issues are relegated to 
appendices. 



2. Electron Stars in the Lowest Landau Level 

Our starting point for constructing magnetised electron stars is the Einstein-Maxwell- 
Dirac action with negative cosmological constant. 



S = J d'x^g [ JL (r+^)~ ^ V - ^ « TaV » 



(2.1) 



Here the field strength is F = dA, while ip is a four-component Dirac fermion 1 . Our 
solutions will be asymptotically AdS 4 with radius L. We choose radial coordinate, r, 
such that the boundary lies at r = 0. 

We impose standard boundary conditions on all fields. For the fermions, this means 

(1 - F> = at r = (2.2) 

while the gauge field has Dirichlet boundary conditions. Our interest lies in dual 
boundary theories at finite density, with chemical potential //, and in the presence of 



1 Our conventions: The covariant derivative contains both spin and gauge connections, "D^ip = 
(dp — iA^ + ^u} lli b c \r a ,T b ])i(). The vierbein e£ translates between tangent space indices, a, b and 
spacetime indices /U, v. Our gamma matrices carry tangent space indices and obey the Clifford algebra 
Clifford algebra {T Q ,r b } = 2ry afc . We define the chiral gamma matrix as T 5 = -ir*r T T a T J '. The 
gamma matrices have the hermiticity property r*(r a ) 1 T t = T a and conjugate spinors are defined by 
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a magnetic field, B. These are implemented through the requirement that the gauge 
fields asymptote to 



A t -¥ —\i , A y ->■ Bx as r -> (2.3) 

Here we have chosen Landau gauge for the magnetic field. As usual, the gauge field is 
not translationally invariant but all gauge invariant quantities, such as the field strength 
F xy , are invariant under translations. In this paper, we restrict our attention to such 
translationally invariant solutions, a seemingly benign fact but one which will have 
important consequences later. The most general form of the metric is 

I 2 

d, s 2 = — (-f(r)dt 2 + h{r)dr 2 + dx 2 + dy 2 ) . (2.4) 

where f(r),h(r) — > 1 as r — > 0. The temporal gauge field also gains a radial dependence, 
A t = A t (r). However, in the absence of any magnetic monopoles in the bulk, the 
Maxwell equations require that the magnetic field is constant 2 . 

The equations of motion arising from (2.1), subject to the boundary conditions (2.3), 
admit a family of well known solutions given by the dyonic AdS Reissner-Nordstrom 
black hole. In this paper we work only at zero temperature. The metric for the extremal 
Reissner-Nordstrom black hole takes the form (2.4) with 



1 4r 3 3r 4 , , 

f r = TT~, = 1 - — + — 2 - 5 
h[r) r h r l 



where the black hole horizon, r^, lies at the root of the quadratic 

6e 2 L 2 



B 2 rt + ^r 2 h — =0 (2.6) 



K 



The temporal component of the gauge field has the profile 

(r 
1 
r h 

Importantly, the fermions are not excited in the Reissner-Nordstrom background. 
Yet this need not be their ground state. As first shown in [6], and explored in some 
detail in [7, 8, 9], for electrically charged black holes the local bulk chemical potential 
may be large enough to excite a bulk fermionic charge density which screens the electric 
field F tr and backreacts on the geometry. The resulting configuration is the electron 
star. The goal of constructing an electron star is to find further solutions f{r), h(r) 
and A t (r) supported by fermions. 



2 Physics in the presence of bulk magnetic monopoles was studied in [16]. 
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As explained in the introduction, a full understanding of the electron star has only 
been achieved in certain regimes of parameter space. This is best developed in the limit 
mL ^> 1, where the Compton wavelength of the fermions is much smaller than the 
curvature scale of spacetime. Here a Thomas-Fermi approximation can be applied for 
the density of states and the fermions are subsequently treated as a perfect, charged 
fluid coupled to gravity in the usual Tolman-Oppenheimer-Volkov manner. Such a 
method was also applied to study neutral (neutron) stars in AdS [17, 18]. 

Electron stars have also been studied in the presence of a magnetic field. Small fields, 
B <C /i 2 — which is relevant if comparing to, say, strange metal phenomenology - were 
studied in [8] where magnetic oscillations, including the Kosevich-Lifshitz formula, were 
recovered as quantum effects in the bulk. More recently, the effect of magnetic fields on 
electron stars in theories with dilaton couplings were discussed in [19]. Both of these 
works treat the electron stars in the fluid, Thomas-Fermi approximation. 

The purpose of this paper is to show that there is another regime where the problem 
of electron stars becomes tractable. This is the regime of large magnetic fields 

B > /i 2 (2.7) 

Here the physics is dominated by the lowest Landau level, both in the bulk and, in an 
appropriate sense, in the boundary. Indeed, one could say that focus of this paper is to 
understand what it means to be in the lowest Landau level of the strongly interacting 
boundary theory. 

Let us first recall some simple facts about free, massless Dirac fermions in d = 2 + 1 
dimensions. With a magnetic field turned on, the spectrum rearranges itself into a 
tower of relativistic Landau levels, with energies E n = V2Bn, for n = 0, 1, 2, . . .. 

The lowest, n — 0, Landau level has vanishing energy and is spin polarized: the 
usual +B/2 zero point energy of the lowest Landau level is cancelled by a —B/2 
Zeeman splitting for the spin down electrons. The net result is that the degeneracy of 
electrons in a plane of area A is BA/2tt for the n = level. All higher Landau levels 
have degeneracy BA/tt. 

For free fermions placed in a magnetic field B > fi 2 /2, the lowest Landau level is 
filled, with all others empty. As we will see in Section 3, the story is much richer for 
the the strongly interacting d = 2 + 1 dimensional theory that lives on the boundary of 
AdS4. In the remainder of this section, we will show how to construct a bulk electron 
star in the regime (2.7). For such large magnetic fields, it is possible to make progress 
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even when the Compton wavelength is comparable to the AdS radius, mL ~ 1. In other 
words, we construct a quantum electron star, albeit restricted to its lowest Landau level. 



2.1 The Landau Levels of the Electron Star 

Our goal is to construct the bulk electron star in the regime ji < yB. To proceed, we 
decompose the bulk fermions into Landau levels in the x — y plane, each of which is still 
free to move in the radial direction. Aspects of fermions in magnetic AdS geometries 
have been previously discussed in [20, 21, 22, 23, 24]. 

With one eye to later bosonization, it will prove useful to first change coordinate 
system. Although the metric (2.4) is more familiar, the 'tortoise' coordinate, f, defined 
such that the (f — t) plane is conformally flat, is calculationally simpler. We therefore 
work with the metric 

ds 2 = Q 2 (r){-dt 2 + df 2 ) + T?(r){dx 2 + dy 2 ) (2.8) 

where 
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n t (f) = f(r)-^ , ^(f) = - (2.9) 



dr y f(r) ' r 2 ' r 

It is a simple matter to write the action (2.1) for the Dirac fermion ip propagating 
in this background, together with a gauge potential A t {f) and bulk magnetic field 
B. However, when dealing with fermions in curved spacetime, it is often simplest to 
perform a conformal transformation, which amounts to a simple rescaling of the fermion 

In these variables, the Dirac action (2.1) takes the form 

Soirac = - J d 4 x # (r r d f + T^dt - iA t ) + | (v x d x + V{d y - iBxf) - mQ^j jt 

Note that the kinetic terms in the f and t directions are those of a fermion in a flat 
space. 

Because Q and E depend on radial position only, the resulting Dirac equation for ip 
is separable: one can solve for the x and y parts of the separable solution and write 
the general solution as a superposition of all such modes. This provides the familiar 
Landau level decomposition of the bulk fermions, 



/it, OO 
— e~ ikv (X n -i(x, k)P + + X n (x, k)P.) $ n , k (r, t). (2.10) 

n=0 
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The matrices P± are projections onto spin-up and spin-down degrees of freedom. 



P± = -(1 ±iT x T y ) (2.11) 

while the profile functions X n (x, k) are the usual, orthonormal wavefunctions of the 
harmonic oscillator expressed in terms of Hermite polynomials H n , 



and X-i = 0. Each n — 0, 1, 2, . . . labels a separate relativistic Landau level. Mean- 
while, the degeneracy within each level is captured by the integral over k. Individual 
states labelled by k are not translationally invariant in both x and y directions. How- 
ever, translationally invariant states can be constructed by taking a superposition of 
all k states. Because, by definition, X_i = 0, the n = Landau level contains only 
spin-down states. This reflects the fact that we mentioned earlier: the lowest Landau 
level is spin polarised and contains only half the states of higher levels. 

It is common practice to replace the integral over degenerate /c-modes with a discrete 
sum. This is achieved this by temporarily restricting the system to a finite transverse 
area A. It is simple to show that the degeneracy of each spin state in a Landau level is 
BA/2tt, and we thus make the replacement 

„ BA/2* 

k 

Within each Landau level, the remaining dynamics then takes place only in the radial, 
f, direction. We emphasise this by introducing two-component spinors, £ + and 
representing individual spin-up and spin-down degrees of freedom respectively, 

£-n,*(r,t) = P-4> n , k (r,t), ^ +n , k (r,t)=zT 5 Typ + ^ k (r,t). (2.13) 

To accompany these, we define two-component gamma matrices 

Y = r r p_, 7 * = r'p_, 7 3 = r 5 p_ (2.14) 

which satisfy the 1 + 1-dimensional Clifford algebra {7^,7^} = 2rf v , with 7 3 = 7*7 r . 
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After this whirlwind of redefinitions, we can finally write our d = 3 + 1 dimensional 
bulk fermion action as a sum over effective d — 1 + 1 dimensional actions for spin-up 
and spin-down fermions in each Landau level. Substituting the solution (2.10) into the 
bulk action yields 



S 



Dirac 



BA/2-k 

E 



drdt 



^2 i£,-n,k(l r df + i\d t - iA t ) - mQ)£- n , k 

n=0 

oo 

+ ^2 i ^.+n,ki.l r df + 7*0% - iA t ) - mVL)S, +n ,k 

n=l 

oo 



n=l 



2BCl 2 n 



(i^+n,k1 3 ^~n,k - i£-n,kl 3 £,+ n,k) 



We stress again that the lowest Landau level, n — 0, contains only spin-down £_ states. 
Notice that, for the lowest Landau level, the only effect of the curved spacetime lies in 
the position dependent mass in the first term, capturing nothing more than the familiar 
gravitational red-shift. The last term in the above expression describes the energy cost 
to excite a mode in a higher Landau level. 

2.2 Bosonization of the Lowest Landau Level 

Our focus in this paper is on dynamics of the bulk fermions at energies E < yB. 
In this regime, all higher n > 1 Landau levels decouple and an effective dimensional 
reduction takes place, with the fermion dynamics well captured by a simple d — 1 + 1 
dimensional action for the lowest Landau level. The effective dimensional reduction of 
the lowest Landau level plays an important role in number of settings including, for 
example, the phenomenon of magnetic catalysis. (See [25] for a recent review). 

Since we focus on the lowest Landau level, from now on we will drop both the 
n = label and spin — label from the associated spinor and write simply = £_o,fc- 
Moreover, because we are interested in translationally invariant solutions, we are at 
liberty to also write the Maxwell action in a similarly dimensionally reduced form. The 
coupled system of fermion and gauge fields is then described by the action 



S — 5*LLL + ^Maxwell 

with the lowest Landau level fermions governed by 

BA/2-k 



(2.15) 



(2.16) 
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and the gauge fields dictated by a two-dimensional Maxwell action in curved space 
which, for reasons that will become apparent later, we normalize as 



^weii = ~f d 2 x {^F^ + ) (2.17) 



In both actions (2.16) and (2.17), \x and v indices run over f and t only and indices are 
raised with the 2D Minkowski metric rf" . Similarly the integral d 2 x = dfdt. Notice 
that, viewed as an action governing the 2d gauge field A^, only the first term is relevant. 
The second term, involving the magnetic field, plays a role only when determining the 
gravitational backreaction. We have introduced the effective two-dimensional coupling 
constant 

1 1 2vrS 2 

(2.18) 



g 2 (r) e 2 B IF 

From (2.9), we have E 2 /f2 2 = l//(r), the black hole emblackening factor. This means 
that g 2 is constant in an AdS spacetime, but approaches zero close to a black hole 
horizon. 

The action (2.15) captures only the interaction between fermions and gauge fields; 
we will discuss the backreaction of these fields on the metric itself shortly. However, 
(2.15) is very familiar: it is the Schwinger model in d = 1 + 1 curved spacetime, albeit 
with a large degeneracy, BA/2ir, of identical species of massive fermions. As observed 
previously in [26], when written in flat space, with conformally rescaled fermions, the 
effect of the spacetime curvature is captured by a position-dependent effective mass 
mfl and a position-dependent effective gauge coupling g 2 . 

Bosonization in Curved Spacetime 

Fermions in two dimensions have a magical property: they can be treated as bosons. In 
flat space, a two-component Dirac fermion operator £ maps into a real scalar operator 
<p (see, for example, [27, 28] for pedagogical introductions). The kinetic terms of each 
field are related by 

= ^0^0 (2.19) 
The vector current of the fermion maps onto a topological current of the scalar 

£A = ( 2 - 2 °) 
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Meanwhile, the axial symmetry of the fermion becomes a shift symmetry of the scalar, 
(f) — > (p + c, with the currents related by 

eWe = ^<j> (2.2i) 

As we will see shortly, this has an interesting consequence for the 2d and, for us, the 
4d anomaly. Finally, the fermion mass term becomes a sine-Gordon potential in the 
scalar theory, 

— 77iA 

im^ = cos0 (2.22) 

which clearly breaks the shift symmetry, as it should. Here A is a regularization scale, 
into which many subtleties of the bosonization map are swept. We shall discuss some of 
these below. This regularization scale is needed in the definition of the massless field 
since the propagator in two dimensions is a logarithm. It also appears in the definition 
of normal ordered products of fields which we have implicitly used in the expressions 
above. As we shall describe in more detail later, when computing physical quantities 
using the semi-classical bosonized action, A should be identified with an appropriate 
physical scale in the problem. 

For our purposes, there are two further issues that arise in the bosonization procedure. 
The first is because we are working in a curved spacetime. This means that all energy 
scales - and this includes the regularization scale A - must be appropriately red-shifted 
as we move in space. For this reason, we should replace A in (2.22) with 3 A — > Af2(f). 

The second issue involves the fact that we have a large number, BA/2ir, of fermions 
that we wish to bosonize. These are labelled by the index k in (2.16). This results in 
an equally large number of scalar fields <pk- However, individually, each of these states 
breaks translational invariance in the x — y plane. Translationally invariant solutions 
only arise if each of these bosons moves as one. For this reason, we identify each 0^ = <fi. 
The net result is that, when restricted to solutions which are translationally invariant 
in the x — y plane, the action (2.16) is equivalent to 

Slll = -— J d 2 x I — W M + (1 - cos0) + (2.23) 



3 This observation seems to have been missed in a number of earlier papers on bosonization in 
curved spacetime. A correct discussion can be found in [29]. Failure to make the regulator position 
dependent in this manner results in a number of pathologies including, as previously pointed out in 
[30], a breakdown of diffeomorphism invariance and, relatedly, the inability to construct a 4d stress 
tensor. Moreover, bosonization in AdS spacetime only respects the scaling symmetry if the red-shift 
factor is correctly treated. (For us, this scaling symmetry is manifest in the equation of motion (3.6) 
for the scalar). 
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Here, e tr = —e rt = +1. As promised, we will say more about the relevant scale A in 
due course. 

The idea that, in certain circumstances - often involving the presence of a magnetic 
field - it is profitable to bosonize fermions in higher dimensions is not new. Perhaps the 
first application - and the one which provided inspiration for the current work - is the 
solution of s-wave scattering of fermions off a magnetic monopole due to Callan [31]. 
Indeed, the s-wave fermi zero mode is morally equivalent to the lowest Landau level 
states that we consider here 4 . Yet more recently, there has been a suggestion to use 2d 
bosonization in conjunction with Eguchi-Kawai reduction to study the planar limit of 
certain 3d field theories [32] . There was also a recent interesting proposal for a genuine 
bosonization of fermions coupled to Chern-Simons terms in d = 2 + 1 dimensions [33] . 

2.3 Fermions, Kinks and Anomalies 

In the following section we will study solutions to (2.23) in asymptotically AdS space- 
times in some detail and explore their implications for the boundary field theory. Here 
we just make a few general comments on the physics captured by this action, all of 
which is standard fare in the bosonization literature. This will also shed some light on 
the meaning of the scale A. 

First, let us ignore the coupling to the gauge field and focus only on the first two 
terms in (2.23). The potential admits vacua at each = 27m, for n e Z. Kinks which 
interpolate between adjacent vacua have unit topological charge and hence, from (2.20), 
are identified with the fermions of the original theory. This identification also tells us 
that — in the absence of any gauge interaction — we should identify A ~ m, so that the 
mass of the fermion matches that of the kink. When viewed from the perspective of our 
original theory in AdS4, the kinks are domain walls of fermionic charge, translationally 
invariant in the x — y direction. They have energy ~ BAmQ/2-K. The interpretation of 
this is clear: the domain walls describe BA/2ir bulk fermions, each of mass m, suitably 
redshifted. We will later identify these as the filled lowest Landau levels of various 
carrier bands in the dual theory. 

Now let's see how this story is affected by the coupling to the gauge field. The 
equation of motion for the gauge field shows us that the kinks are electrically charged 

4 This analogy can be seen by starting with a magnetically charged black hole in global AdS. Upon 
taking the limit to the planar Reissner-Nordstrom in the Poincarc patch, the s-wave zero modes evolve 
to the lowest Landau level. 



12 



as expected, 




1 



(2.24) 



This is where issues of the regularization scale A become subtle since, in some contexts, 
it is necessary to take A proportional to the gauge coupling g when treating (2.23) semi- 
classically. In the next two paragraphs we will explain why we shouldn't do this and 
must instead retain A ~ m. Readers uninterested in these subtleties are invited to 
jump to equation (2.25). 

Perhaps the simplest, and most compelling, way to determine the regularization scale 
A is to use the fact that we're working in a higher dimensional, curved spacetime. In 
Appendix A, we construct the fermionic contribution to the d — 3 + 1 stress tensor 
in terms of the bosonized field 0. This stress tensor depends explicitly on A and the 
requirement that is covariantly conserved puts strong constraints on what A can be. 
This is essentially because the effective gauge coupling g 2 , defined in (2.18), depends 
on the radial position on a way which remembers its four- dimensional origin. We show 
that the only consistent choice of regularization scale is A ~ m. 

It is also useful to address the issue directly in a two-dimensional language. In the 
usual discussion of the Schwinger model in flat space, one puts A ~ g rather than 
A ~ m. (At least this is true in the regime g 2 ^> m 2 where the bosonized action is 
weakly coupled). Here the choice of A reflects the fact that the mass of the electron 
is renormalized by the interaction with the gauge field, as can be seen in the bosonic 
picture by integrating out the gauge field. However, even here one should only put 
A ~ g for the Schwinger model with a single species of fermion. In the presence of two 
fermion species, £i and the story changes. While the singlet boson <pi + <p 2 acquires 
an effective mass proportional to g, the other boson combination 0i — 02 acquires a mass 
proportional to w?l z g x l z . A careful treatment shows that the correct normal ordering 
scale is now A ~ m 1 ^ 3 g 2 ^ 3 [34]. Extending this argument to Nf species of fermions, the 
normal ordering scale becomes A N f +1 ~ m N f~ 1 g 2 [35]. Since we have BA/2ir species 
of fermions, we should take the Nf — > 00 limit. We again find A ~ m as promised 5 . 

5 There is yet a third way to fix the regularization scale A which, moreover, provides insight on a 
slightly different topic. One physical quantity in which the regularization scale appears is the chiral 
condensate = A/7r and, in the lowest Landau level, it is known that the chiral condensate is 

proportional to the mass (see, for example, [36]). Furthermore, with the identification A ~ m, it is 
simple to show that a number of higher dimensional phenomena which rely on the lowest Landau level 
- such as the chiral magnetic effect and the chiral spiral (see again [36] for a review) - follow trivially 
from the equations of motion of the bosonized theory. 
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The arguments above for A ~ m do not fix an overall, order one coefficient. For this 
reason, we'll replace mA — > m 2 R , the renormalized mass, in the above expressions. The 
equation of motion for then reads 



The second term on the right-hand side is particularly interesting. It captures the axial 
anomaly in two dimensions, as can be seen by substituting the expression (2.21) for the 
axial current. However, our theory has a four dimensional origin and all expressions 
can be rewritten in in terms of the four dimensional fermion ip. Setting m = in (2.25) 
and taking into account the degeneracy, BA/2ir, of the lowest Landau level, we see 
that the two-dimensional anomaly is, at heart, the four dimensional anomaly 6 



where the axial charge is Qaxiai = ijjT t r 5 ip. The axial anomaly is a property of the 
quantum theory of the fermion field. It is a beautiful feature of bosonization that the 
classical equations of motion capture this quantum aspect of fermions. 

Gravitational Backreaction 

So far, we have restricted our discussion to the dynamics of fermions and gauge fields in 
the background of a fixed, albeit arbitrary, asymptotically AdS4 background. Ideally, 
one would like to generalise this to include backreaction on the metric itself and, indeed, 
in the framework of bosonization there is neither conceptual nor technical obstacle 
to this. In appendix A, we compute the four-dimensional stress tensor arising from 
fermions in the lowest Landau level and reduce the Einstein-Maxwell-Dirac equations 
to a set of coupled, ordinary differential equations. 

Nonetheless, in this paper we do not solve the full backreacted equations. The reason 
for this is rather trivial: in the regime, B > /z 2 in which we are working, the backreaction 
on the metric is dominated by the magnetic field B. Furthermore, while the electron 
star screens the electric field in the bulk, it does nothing to change the magnetic field. 
The upshot of this is that, at least in the regime where the most interesting things 
happen, the fully back-reacted solutions do not differ substantially from the purely 

6 Hcrc wc include only the contribution from the lowest Landau level. Higher Landau levels trans- 
form as £„± — > e^ ldl / 2 £± under the chiral rotation ip —> e l6T ^ 2 ip. The axial anomalies from spin-up 
and spin-down states of the higher Landau levels cancel out, leaving only the contribution from the 
lowest Landau level. 



d^d^ = m 2 R n 2 sin + e^F^ 



(2.25) 
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magnetic Reissner-Nordstrom black hole background, described by (2.5) and (2.6) with 
p = 0. For this reason, in the following section we use bosonization to study the lowest 
Landau level of electron stars in the purely magnetic Reissner-Nordstrom black hole. 

3. Phenomenology of the Lowest Landau Levels 

In this section we describe in some detail the solutions to our bosonized system and their 
implications for the dual boundary field theory. We will focus entirely on equilibrium 
properties but will, nonetheless, find a rich phenomenology. 

One of the main obsessions in this section is on the charge density, p, of the boundary 
field theory. As always, this has the interpretation of the bulk electric field near the 
boundary. More precisely, working in A r = gauge, the charge density can be read off 
from boundary behaviour of A t , 

A t -> -p + e 2 pr as r -»■ (3.1) 

From this we see that the charge density p is proportional to the value of the electric 
field F tr on the boundary. This charge density can have two qualitatively different 
sources [37, 38]. Firstly, the electric field can be sourced by charged matter in the bulk. 
For us, this charged matter means the fermions or, after bosonization, the kinks of 
the scalar field. In the boundary theory, this corresponds to charge carried by gauge 
invariant fermionic operators ^ dual to the bulk field ip. Such charge is referred to 
as 'mesonic' or 'cohesive'. As shown in [13], this charge contributes to the Luttinger 
count. 

An alternative source of electric field is provided by the horizon of a Reissner- 
Nordstrom black hole. In the boundary theory, such charge is carried by the elementary, 
non gauge-invariant fields of the large N gauge theory (whatever they may be). Such 
charge is said to be 'fractionalised' 7 . 

Much of the focus in the story of electron stars is to determine which of the two 
sources of charge - cohesive or fractionalised - is the stable ground state in any given 
setting. Indeed, one can find situations where both play a role [39]. The main purpose 
of this section is to describe how these two types of charges arise in the presence of a 

7 A particle theorist would call such charge dcconfined since particle theorists are used to having 
fundamental non-Abelian gauge fields at a high scale. The term fractionalised is more appropriate 
in the condensed matter setting where non-Abelian gauge fields first emerge at some low scale then 
de-confine in a non-trivial manner, the net effect being that elementary particles - such as the electron 
- appear to fractionalise. 
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strong magnetic field. We will find that there is an interesting interplay between the 
two. 

Although we are ultimately interested in the physics in the background of a bulk 
Reissner-Nordstrom black hole, we will first start in a simpler setting where we can 
better illustrate the main ideas of formation of an electron star in the framework of 
bosonization. In Section 3.2, we will place our electron stars in a black hole spacetime 
and illustrate a number of different phenomena that occur in the presence of a horizon. 

3.1 Electron Stars in a Hard wall 

To describe the formation of an electron star, we start with the simple situation of a 
AdS hard wall. The metric is 

L 2 

ds 2 = —(-dt 2 + dr 2 + dx 2 + dy 2 ) (3.2) 

The UV boundary lies at r = 0. In the infra-red, we assume that the space stops 
abruptly at r = r*. This is the the hard wall. This geometry was also used in [13] in 
the construction of a quantum electron star. 

We thread this geometry with a constant magnetic field F xy = B. For the purposes 
of this toy model, we will assume that this magnetic field is emitted by the hard-wall. 
Furthermore, we ignore the back-reaction of the magnetic field on the metric. (Such 
ignorance is justified if Br\ <C ^Je 2 L 2 / k 2 ). 

Already in vacuum, there are a number of interesting effects experienced by fermions 
in a the magnetic hard wall geometry [23, 24]. Here we differ from these papers by 
turning on a charge density in the boundary field theory. This means that, close to the 
boundary, the bulk gauge field takes the form (3.1) with p ^ 0. We will further insist 
that no electric field is emitted from the hard-wall; the electric field must therefore be 
entirely sourced by bulk fermions which, after bosonization, means kinks of the scalar 
field <j). We would like to ask what charge densities p can be built in this manner. 

In fact, this is a variant of a classic problem in physics. If you stand on your head, 
the kinks which source the electric field can equally well be thought of as screening 
an electric field emitted from the boundary. This is the question of charge screening 
in the Schwinger model. Famously, after bosonization this classic problem becomes a 
classical problem 8 [15]. 

8 This statement is almost true. The bosonizcd theory can be treated semi-classically in the regime 
g 2 3> m 2 R Q 2 where its self-interactions are weak. This translates to the requirement Be 2 3> m 2 L 2 /r 2 . 
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Working in Coulomb gauge, and assuming a time independent configuration, the 
classical equations of motion of the scalar and gauge field are 

e 2 B m 2 L 2 

d 2 A t = -d r <p and d 2 <p = — ^— sin0 - 2d r A t (3.3) 

4n 2 r z 

The first of these equations, Gauss' law, can be integrated trivially to give an expression 

for the electric field in terms of the scalar field. 

3 r A t = -0(0-0 o ) (3-4) 

Here we have introduced an integration constant, 0o- To understand the meaning of 
this integration constant, we need to discuss boundary conditions for <f) at r = 0. These 
are inherited from the boundary conditions (2.2) for the fermions 9 and are simply = 
[31]. Comparing (3.4) to (3.1), we learn than 0o determines the charge density 

p = — 2 h (3.5) 

In the context of free electrons, </>o/27r is called the filling fraction] it tells you how 
many Landau levels are filled. In our case, we are already restricting ourselves to the 
lowest Landau level but, of course, we are describing a strongly interacting, large N 
theory. We will continue to refer to </>o/27r as the filling fraction. However, it will now 
have a slightly different interpretation, telling us how many bands - or different species 
- of fermions have their lowest Landau level filled 10 . 

Substituting (3.4) into the second equation in (3.3), we arrive at a second order 
differential equation for the scalar, 

9 r 2 0=^sin0+ J-(0-0o) (3.6) 



s 

where we have introduced a new length scale, 

r s = (3-7) 



e 2 5' 

This represents the characteristic length scale over which the bulk electric field is 
screened (at least in the case of small mass). 

Notice that, for m/0, a semi-classical treatment of the bosonized action is not strictly valid close to 
the boundary. As we shall see, at least for small masses the most interesting physics occurs suitably 
far from the boundary. 

9 In fact, there appear to be interesting choices in the allowed boundary conditions for both fermions 
and scalars. We have relegated the details to Appendix B. 

10 A comment on the regime of validity: we cannot increase <j>o too much or we will start to fill higher 
Landau levels and our framework of bosonization will break down. We are parametcrically safe, with 
[i < yB, provided that </>o < 1/e. This means that we can fit up to n ~ 1/e kinks in our geometry. 
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Before solving (3.6), we must first discuss the boundary conditions which we impose 
at the hard wall. The canonical choice is to again impose (1 — T r )ifr = at r = r*, just 
as we did at the boundary (2.2). In the UV, this translated into the requirement = 
and one might think that that we should also insist on this at the hard wall. But this 
is not quite right. The boundary condition (2.2) is invariant under a chiral rotation of 
2tt. Yet, under such a rotation, cfr — > cfr + 2ir This means that the boundary condition 
on the fermion actually allows a discrete choice of boundary conditions on the scalar, 
namely 

(f) = 27m, n G Z at r = r* (3-8) 

The one further requirement that we need to impose is that the hard wall emits no 
electric field. In the bosonized language, the relationship (3.4) tells us that this becomes 
cfr = <fro at r = ?v Combined with (3.8), we learn that the hard wall background admits 
translationally invariant electron stars in the lowest Landau level only for integer valued 
charge density 

B 

The same result was found in [40] which examined the Luttinger count for electron 
stars in a hard wall and magnetic field. If one tries to impose a non-integer filling 
fraction, the resulting state is presumably not translationally invariant in this setting. 
We will see shortly that this situation changes in the presence of a black hole horizon. 
(In Appendix B we also describe boundary conditions at the hard wall which result in 
fractional charge density). 

For massless bulk fields, m = 0, the solution to (3.6) is a simple exponential 

1 — e~ r l Ts 

4> = i _ e _^ lTs (3.9) 

This describes the profile of fermion charge density in the bulk sourcing the electric 
field. 

For m ^ 0, we solve (3.6) numerically 11 . The resulting profiles are shown in Figure 1. 

Since the mass makes it unfavourable to place fermions near the UV, the charge density 

of the electron star is pushed further away from the boundary as the mass increases. 

Most strikingly, the electron star develops kinks of localised charge density within 

its profile. These are simply the kinks of the sine-Gordon model, now transplanted 

into AdS. If you allow more fermions in each Landau level by increasing B, then the 

boundary electric field increases and the kinks are drawn towards the UV boundary. 

n We use a fourth-order Runge-Kutta algorithm with variable step size. The IR boundary condition 
is implemented using a shooting method. 
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Figure 1: Electron stars in a hard wall. The solid line shows the bulk fermionic charge 
density, normalised as Jt = r s d r (j)/2TT as a function of radial position f = r/r+; the dotted 
line shows the bulk electric field, again normalised as E = r 2 F r t/i:. The left hand plot has 
filling fraction (fio/2ir = 1 with just a single kink; the right-hand plot has 0o/2tt = 3 and three 
distinct peaks. Both of these plots are made with m 2 R L 2 = 16 and r 2 /r 2 = 1/10. 

Within each kink, <fi increases by 2ir and, using (3.3), therefore sources a charge 
density in the boundary of B/2n. This is the charge density obtained by filling a 
lowest Landau level of states. From the perspective of the boundary field theory, 
we are filling the lowest Landau levels of different species - perhaps a better name is 
different carrier bands - of fermions. In the bulk, these bands are associated to different 
radial harmonics of the fermion field 12 . Note, however, that unlike in the free theory, 
the Landau levels are not delta-functions in energy space. They have been broadened 
by the strong interactions of the boundary theory, captured by the AdS geometry. 

Viewing the radial coordinate r in AdS as an energy scale, E ~ 1/r, it is tempting 
to view the kink profiles in Figure 1 (and those of Figure 2) directly as a plot of the 
density of states of the lowest Landau level of the boundary field theory. 

3.2 Black Holes 

We now turn to the main case of interest: electron stars in the background of a magnetic 

12 The presence of different carrier bands is also seen in the original electron star papers where it 
results in concentric rings of Fermi surfaces [9, 12]. Indeed, in the limit mL 3> 1 necessary to treat 
the electron star as a fluid, these Fermi surfaces are densely packed. The same phenomenon occurs 
here: the width of a kink at position r is roughly r/m. As m increases, the kinks become sharper and 
more closely spaced. 



19 



Reissner-Nordstrom black hole. We take the geometry to be 

L 2 ( dr 2 
ds 2 = — -f{r)dt 2 + -r-- + dx 2 + dy 2 
r 2 V f( r ) 



with 



4r 3 3r 4 



3 1 4 



The position of the horizon is fixed by the magnetic field, 



6e 2 L 2 1 



k 2 B 

Our use of this background requires some explanation. As we shall see, we will be 
interesting in solutions containing both cohesive charge, carried by the star, and frac- 
tionalised charge lying behind the horizon. Yet we do not allow the position of the 
horizon to depend on the electric charge of the black hole. This sounds inconsistent. In 
fact, as we show in Appendix A, the backreaction of the electric field on the position of 
the horizon is negligible provided e 2 0o *C 1 and e 2 4>omRL <1. It is easy to understand 
why. Even in the UV, the magnetic field B is larger than the electric field. By the time 
we get down to the horizon the vast majority of the electric field has been screened by 
the star and what remains does not meaningfully affect the position of the horizon. 

The equations of motion again tie the bulk electric field to the scalar, 

«9 r A t = -0(0-0 o ) (3.10) 
Meanwhile, the second order equation governing the scalar itself is now given by 



d r (f(r)d r <f>) = ^^sin0 + ^(0-0o) (3.11) 



where, as in the previous section, the screening length is defined as r 2 = 27r 2 /e 2 B. 

The primary difference with the hard wall background lies in the boundary conditions 
imposed on 0. In the UV, we again set = at r = 0. As before, this allows us to relate 
0o to the charge density: p = B^o/An 2 . However, we do not impose Dirichlet boundary 
conditions in the IR. Instead we require only that (ft and <9 r are regular on the horizon. 
However, because f(r) has a double zero in the extremal Reissner-Nordstrom black 
hole, this in turn requires that <j){rh) lies in an extremum of the effective potential 13 

2 

V IR {<j>) = m\L 2 {\ - cos0) + ^(0 - O ) 2 . (3.12) 
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Figure 2: Typical electron star configurations in the black hole. The solid lines show the 
charge density, Jt = fr s d r (f)/2Tr as a function of radial position r = r/r^. The dashed lines 
show the normalised electric field, E = r 2 F r t/ir. The left-hand graph has filling fraction 
4>o/2ir = 2 and is plotted rn 2 R L 2 = 4 and r\jr 2 = 5. Note that only a single peak is visible; 
this is typical for low values of the mass. The right-hand graph is plotted with m? R L 2 = 16 
with r^/r 2 = 10 with a fractional filling number </>o/27r = 3.25. Note that this results in a 
residual electric field at the horizon of the black hole. 

The first term is due to the mass of the fermion. The second term arises after integrating 
out the gauge field and can be thought of as capturing the electrostatic energy between 
fermions. 

We again solve (3.11) numerically. Results for typical values of the parameters 14 
are plotted in Figure 2: solid lines denote the local bulk charge density, dashed lines 
the electric field. These plots look very similar to those in Figure 1 because most 
of the kinking is taking place in the AdS4 region of the geometry. However, there 
are differences. Most importantly, the allowed charge density is no longer restricted to 
integer filling fraction. Solutions with fractionally filled Landau levels happily exist and 
we have presented one in the right-hand graph. Such solutions are typically (but not 
always) accompanied by a residual electric field at the horizon, showing the existence of 
fractionalised charge. We will devote Section 3.4 to a deeper analysis of the relationship 
between fractional filling and fractionalised charge. (Just because they're both called 
named after fractions doesn't mean they're the same thing!) 

13 A full analysis of the requirements in the near horizon AdS2 x R 2 regime is given in Appendix C. 

14 These plots are made in the regime r s -C or, equivalcntly, e?L 3> n. This ensures that the 
screening, and hence kink formation, takes place before we hit the horizon. In the opposite regime, 
r s 3> r/j, the electric field is largely unscreened for much of its profile, although some screening docs 
take place near the horizon. 
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Figure 3: A graph of how the normalised charge density in the boundary theory, p = pr^/ir = 
4>o/2ir, varies with the chemical potential, p = pr s /ir. Note that the vertical axis corresponds 
to filling fraction. The normalization is such that in this, and all other plots, the massless 
fermion gives rise to a straight line of unit slope. This graph was made with m 2 R L 2 = 6 and 



3.3 Charge Jumps 

As we will now see, the presence of the kinks in the bulk shows up in a number 
of physical properties of the boundary field theory. Perhaps the simplest physical 
observable is the boundary charge density, p. This can be plotted as a function of the 
chemical potential p, given by integrating the electric field between the boundary and 
horizon 



The relationship between p and p in our system is somewhat intricate and changes 
qualitatively as various parameters are varied. In this section and the next, we take 
some time to describe this. 

For massless fermions, the relationship between the charge density and chemical 
potential is linear. As we turn on the mass m of the fermion and, correspondingly, the 
sine-Gordon potential for the bosonized scalar, kinks form in the bulk. Their effect is 
to induce wiggles in the relationship between p and p. These wiggles occur close to 
integer filling fractions and are largest for small p, dying out as the charge density is 
increased. A typical example is shown in Figure 3. 

It is worth stressing the obvious: the relationship between p and p in this regime is 
continuous. We can contrast this to the situations for free fermions in a magnetic field 



rl/r 2 s = 8. 




(3.13) 
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Figure 4: The total charge density in the boundary theory, p = pr 2 /ir = 0o/27T. This graph 
was made with parameters m 2 R L 2 = 12 and r 2 /r 2 = 20. Note the discontinuous jumps in the 
first five bands. 

which would result in step-like behaviour every time the chemical potential crosses the 
energy of each Landau level. As we have already noted, in our strongly interacting 
theory, the energy levels of the bands are broadened and the steps are smoothened out, 
with only ripples surviving. 

However, as the mass rn 2 R L 2 is increased further, the wiggles in the pvs/i plot become 
more pronounced. Eventually, they grow so large that the graph p(p) is no longer single 
valued. A characteristic example 15 is shown in Figure 4. For values of p which allow 
solutions with two different charge densities p, one must compute the free energies to 
determine which wins. We will discuss the free energy of the electron stars in Section 
3.5 -see equation (3.17) - in the context of quantum oscillations. In the present case, 
the result is simple: the solution with the smaller free energy is that with the smaller 
value of p. This results in a discontinuous jump in the charge density as the chemical 
potential is increased. (A similar quantum Hall-like step was seen for lowest Landau 
level in the D3-D7 probe brane system [41, 42]). 

The discontinuities in the total charge density p occur only for the lowest bands. 
As the mass is increased (all other parameters remaining fixed), more and more bands 
suffer the discontinuity. 

3.4 Fractionalised and Cohesive Charge 

Above, we have seen that there are discontinuities in the total charge density only if the 

This graph is created by working backwards, fixing p and then determining /i. The function p,(p) 
is single- valued, but is not monotonic. 
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mass is suitably large. However, even in the case that p appears smooth - as depicted 
in Figure 3 - if we look more closely we can see that the familiar Landau level steps are 
not as far away as one might imagine. As we reviewed at the start of this section, charge 
density in the holographic framework comes in two forms: cohesive charge, carried by 
bulk fields, and fractionalised charge, hidden behind the horizon. Let us see how the 
charge density p is split between these two options. 

The amount of fractionalised charge is determined by the electric field at the horizon. 
Using (3.10), we can write this as 

Pfrac = (00 ~ <f>(r h )) (3-14) 

From the discussion in the previous section, we know that 4>(rh) necessarily lies at an 
extremum of the effective potential (3.12). 

It is simple to extract the amount of fractionalized charge pf rac from our numerical 
solutions. We find that as the chemical potential is increased, the fractionalised charge 
density oscillates about zero. (Examples can be seen in Figures 5 and 6). Using the 
considerations above, it is straightforward to check that in all cases the amplitude of 
oscillation is bounded by 

Kac| < ^ (3.15) 

Notice in particular that, in the case of massless fermions, all the charge is carried by 
the electron star; none by the horizon. 

To understand the fractionalised charge density in more detail, we need to look at 
the extrema of the potential (3.12). These depend on the relative size of fn 2 R L 2 vs 
r 2 /r 2 . As we now show, this gives rise to a qualitative difference in the behvaiour of 
the charge density. 

When pf rac Changes Continuously... 

We start by considering m 2 R L 2 < r 2 Jr 2 . In this regime, Vm(<t>) has a unique extremum 
for any value of filling fraction <p /27r. This extremum is a global minimum and regular- 
ity at the boundary means that 4>(rh) must sit in this minimum. With no ambiguity in 
the solution, the fractionalised charge is a continuous function of \i. The fractionalised 
charge that contributes to the total charge density shown in Figure 3 is shown in the 
left-hand graph of Figure 5. Note that the fractionalised charge can be both positive 
and negative. 
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Of course, if the fractionalised charge oscillates and the total charge is monotonically 
increasing, the remainder must be made up in cohesive charge. This is shown in the 
right-hand graph in Figure 5. 

Notice that the fractionalised charge density vanishes for certain values of the chem- 
ical potential. Translating this into filling fraction, we find that pf rac = whenever a 
Landau level is fully filled or half filled. This can be understood by looking again at 
the effective potential (3.12). The charge at the horizon vanishes whenever = 0o is 
the minimum of this potential. This holds for = nir with nfZ. Integer filing corre- 
sponds to even n; half-integer to odd n. At half-integer filling, sits at a maximum of 
the (1 — cos0) potential, giving rise to an extra half kink in the bulk. But, by virtue 
of the condition vn? R I? < r^/r 2 , this maximum of (1 — cos0) becomes the minimum of 
the full effective potential (3.12) which includes the quadratic electrostatic term. 

The fact that pf rac = at integer filling is unsurprising. The fact that pf rac = also 
at half-integer filling in this regime is surprising. We do not understand this from the 
perspective of the boundary theory. 

When p frac jumps... 

As we increase the mass past the critical value, m 2 R L 2 > r^/r 2 , things get more inter- 
esting. In this regime, there exist some values of filling fraction 0o/27r for which there 
are multiple extrema of V/\r(0). A full study of the solutions is a little complicated and 
is described in Appendix C. Despite multiple minima, only one value of <j)(rh) gives 
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Figure 6: When ni 2 R L 2 > r 2 /r 2 , both the fractionalised charge density (on the left) and 
the cohesive charge density (on the right) are show discontinuities. Here we have taken 
m 2 R L 2 = 12 and r^/r 2 = 4. When combined the two discontinuities cancel out to form the 
smooth total charge density that does not look qualitatively different from Figure 3. 

rise to a stable field configuration at any filling fraction. This preferred 4>(rh) is al- 
ways a local minimum of Vm((j)), but not necessarily a global minimum. Furthermore, 
as we vary the chemical potential, there are values at which the preferred minimum 
jumps discontinuously. At these points, the fractionalised charge also jumps. A typical 
example is shown in the left-hand graph of Figure 6. 

The discontinuity also arises in the cohesive charge density. A plot is shown in the 
right-hand side of Figure 6 and exhibits the characteristic step features expected of 
Landau levels. (Although these are really different bands, all of which lie in the lowest 
Landau level). These steps continue indefinitely for higher bands. Note, however, that 
the plateaux are not precisely horizontal; the Landau levels in the strongly interacting 
boundary theory are not precisely flat bands. 

We stress that the step-like behaviour seen in the cohesive charge is not obviously 
related to the discontinuities appearing in the total charge that we exhibited in the 
previous section. In particular, the jumps in the fractionalised and cohesive charge 
shown in Figure 6 occur both in the regime where the total charge is monotonic (as 
shown in Figure 3) and also in the regime where the total charge itself exhibits a finite 
number of jumps (as shown in Figure 4). 

Of course, for free fermions the lowest Landau levels are the domain of quantum 
Hall physics. In this paper, we have restricted attention to equilibrium properties of 
the system and have not attempted to compute transport. Needless to say, it would 
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Figure 7: The dimensionless magnetic susceptibility, x = (6e 2 L 2 /5 2 /K 2 ) 1//4 x plotted against 
the filling fraction, 0o/27r. The parameters chosen for this plot are rn 2 R L 2 = 12, r^/r 2 = 20. 

be of great interest to do so, although the techniques of bosonization do not obviously 
lend themselves to the task. Nonetheless, plots such as Figure 6 show that one will find 
oscillatory or step-like behaviour for any process in which the fractionalised and cohesive 
charges contribute differently. In the next section, we present a simple example. 

3.5 Quantum Oscillations 

In this section, we wish to compute how the magnetic susceptibility varies as a function 
of magnetic field B at fixed charge density p. As we increase B, the occupancy of each 
Landau level increases and, correspondingly, the filling fraction decreases.. Each time 
the number of filled bands changes, the magnetic susceptibility oscillates. These are 
the famous de Haas van Alphen oscillations, although in our case we are filling the 
lowest Landau levels of different bands as opposed to higher Landau levels of the same 
band. From (3.5), the period of oscillation is 



For three-dimensional metals, this is often written as A(l/B) = litejAp where Ap is 
the extremum area of the Fermi surface. 

To compute the magnetic susceptibility, we first need to evaluate the free energy 
of the bulk fermions and electric field. Because we are working at zero temperature, 
this simply means the energy. Integrating out the electric field from (2.23) yields an 
effective Hamiltonian for 0, 




(3.16) 




BA 
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Note that this is closely related to the effective potential defined in (3.12). The magnetic 
susceptibility is then given by evaluating the on-shell Hamiltonian and taking the second 
derivative, 

1 d 2 H eS 
X ~ A dB 2 

The result is plotted in Figure 7 and clearly exhibits quantum oscillations. This particu- 
lar plot was made in a regime where both the total charge density and the fractionalised 
charge density vary smoothly (as in Figures 3 and 5). However, the figure looks qual- 
itatively the same for the other regimes. Notice, in particular, that we did not need 
to turn on a temperature to smooth out the oscillations; quantum effects have already 
achieved this for us 16 . 

The free energy computed above does not contain the contribution from the gravita- 
tional background. This contains an overall contribution from the Reissner-Nordstrom 
black hole which, at zero temperature scales as x ~ 1/ V~B, providing an overall enve- 
lope for the contribution above. There will furthermore be an oscillatory component 
due the horizon moving backwards and forwards as the fractionalised charge waxes and 
wanes. However, this contribution depends on a new parameter, k, and cannot remove 
the oscillations from the electron star computed above. 

Quantum oscillations were seen previously in the probe approximation [43, 44], where 
they are a sub-leading 1/N effect, and also in electron stars computed in the fluid 
approximation [8], where the standard Kosevich-Lifshitz formula for their amplitude 
was reproduced. It would be interesting to understand the amplitude and temperature 
dependence of the oscillations in the present case. While we have not made analytic 
progress with this question, numerics suggest that the increase in the amplitude of 
oscillations is fairly well approximated by the \ ~ 1/B 2 behaviour of the Kosevich- 
Lifshitz formula. 
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A. Appendix: Einstein's Equations and the Stress Tensor 

The purpose of this appendix is to construct the 4d stress tensor in bosonized form. 
As we shall see, the requirement that the stress tensor is conserved puts severe con- 
straints on the method of bosonization in curved space and, in particular, the choice 
of regularization scale A. 

With the stress tensor in hand, we will write the Einstein equations for translationally 
invariant states in bosonized form. Although we do not solve these equations, we will 
check that in the regime of interest the gravitational solution is dominated by the 
magnetic field. 

The Einstein equations resulting from the action (2.1) take the form 



1 3 

Rpv ~^RQpv 9 'fif K 



In the conformal (r,t,x,y) coordinates, with metric ansatz (2.8), the various compo- 
nents of the Einstein tensor are 

1 3 3f2 2 
Rtt - xRgu ~ Tl9tt = -X- 2 (d f £) 2 - 2E~ 1 d^E + 2ir 1 £- 1 ,9 r m~£ + — 

2 Lr L z 

1 3 3f2 2 

Rfr — T^Rdff — ~T^9ff = (<9f£) 4~ 2f2 S df-VldfYj — — —— 

2 L 2 L l 

1 3 3E 2 

R xx - ^Rg xx - -^g xx = ^r 2 x«9 2 £ + ir 3 s 2 <9 2 ft - ir 4 £ 2 («9 r -ft) 2 - — 

The stress tensor contains contributions from both the gauge fields and fermions 

rp _ rpgauge I T-rfermion 

~ 1 pV ' 1 pi/ 

The contribution from the Maxwell term can be derived by varying the action with 
respect to the metric 

Meanwhile, the stress tensor for fermions in curved space can be derived by varying 
the action with respect to the vierbein. (See, for example, [45]). The result is 

-i' 



rpgauge _ / p pp p ppa\ 

pv o ^ w L u a p u ) 



t 1 ™ = - A ^ap ( d v + -u uM [r\ r c ] - iA v ) r> + h.c. ) + (im*>v) 
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Since we are interested only in states that enjoy translational invariance in the x — y 
plane, we can write 

1 



T-[a> J dxdy T^ u , 

Substituting in the mode expansion (2.10), and performing the integrals over x and y, 
we arrive at expressions for the stress tensor in terms of the Landau level fermions 

T-ifermion _ 



T-ifermion 



T-ifermion _ 



-y / OO OO ^ 

j^a Yl [Y^-n^idt - iA t )C- n ,k + ^2^+ n ,kl t (d t - iA t )£, +n>k + h.c. 

k \n=0 n=l / 

/ o \ 

k \n=0 n=l / 

1 ( °° 

^2 Yl ( 53(^-n.fcT r (^ ~ iA)£-n,k ~ i^-n,kl td fi-n,k) 
k \n=0 

oo 

+ ^2(i£, +n , k Y(d t - iA t )£ +n ,k - iC+n^dfC+^k) + h.c. 



71=1 



As is clear from these expressions, each of these components receives contributions 
from all Landau levels. However, when we come to compute the T^ rmion and T^ rmion 
components of the stress tensor, only the higher Landau levels contribute, 



T-ifcrmion T-ifcrmion _ 



^ ( i^B^L+l Z in- + h.c. J 

k \n=l J 



This is another manifestation of the fact that fermions in the lowest Landau level 
undergo an effective dimensional reduction and are, therefore, susceptible to bosoniza- 
tion. Since we will excite only fermions in the lowest Landau level, we are interested 
in solutions with T^ rmion = j^ rmion = o. Note that this automatically means that we 
are working in a regime that is far from the isotropic stress energy tensor assumed in 
Thomas-Fermi approximations to the electron star. 

Restricting now to fermions in the lowest Landau level, we would like to construct 
expressions for T| rmion , T~ rmion and T~f mion in terms of the bosonized field <fi. In fact, 
requiring a consistent stress tensor provides very non-trivial constraints on the choice of 
regularization scale A. Specifically, we will ask that our stress tensor is gauge invariant, 
symmetric and, most importantly, conserved when the equations of motion are obeyed 

= 
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These criteria are only satisfied if we pick A ~ m. The resulting expression for the 
stress tensor for translationally invariant and time invariant states is 

1 (B \ ( 1 ,„ ,, 2 m|L 2 S! 2 , \ 1 / 1 ,„ „ l2 , n 2 



T " ~ E 2 \2tt ) WW + ^to-V + 2? l^'' + * B 



£ 2 V27ry V8tt v ^ 4tt v T! ) 2e 2 V^ 2 S 4 

Furthermore, the off-diagonal component vanishes for time-independent states, T ft = 0, 
while, as described above, the pressure in the x and y directions gets contributions only 
from the gauge field, 

1 (Y? 1 
T xx = T yy = — i—(dfA t ) 2 + —B 2 

As a further check of this result, note that when the r — t plane is flat, the stress tensor 
should coincide with the Noether currents arising arising from translational invariance: 
it does. 

Notice that altogether we have five equations of motion, but only four degrees of 
freedom. As usual, there is no inconsistency because the Bianchi identity and the 
conservation of the stress tensor imply that one of the equations is redundant. 

Einstein Equations 

Using the bosonized stress tensor, we can now derive the coupled Einstein-Maxwell- 
Dirac equations in bosonized form, capturing the gravitational backreaction of fermions 
in the lowest Landau level. The resulting equations are simplest if we revert to the r 
coordinates of (2.4) and parameterize the metric as 

ds 2 = % (-e 2C ^D(r)dt 2 + -^-Jr 2 + dx 2 + dy 2 ] 
r i y D{r) J 

Gauss' law can be integrated trivially to give 

e 2 Be c 

d r A t = ^-(<f>-<f> ), 

and the other fields obey a simple coupled system of ordinary differential equations. 
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Background Check 

We now use the results above to determine when we can neglect the backreaction of the 
fermions and electric field on the metric. This holds when the dominant contribution 
to T u is provided by the magnetic field. The relevant equation is 

The phenomena described in Section 3.2 are most clearly observed when r s <C r^. In 
this case, the envelope of is well approximated by the exponential profile (3.9). From 
this, we estimate 

^~ 2 (4> - 0o) ~ \/^ d r<f> ~ ^0 O exp(-r/r s ) 

The kinetic and electric terms in T tt are therefore negligible compared to the magnetic 
term if 

e 2 o < 1 

We must also consider the energy contained within the kinks. Since the kinks are 
embedded in an approximately exponential profile, we estimate that the n th kink is 
located at position rkink given by 

0o - 27m ~ o exp(-r kink /r s ) 

Demanding that the energy of a kink is small compared to the energy of the magnetic 
field, we obtain the condition 

m R Le 2 (f) Q 1 

B. Appendix: Boundary Conditions and Filling Fractions 

In this appendix, we make a few comments about the boundary conditions after 
bosonization. We start in the UV, where we imposed = 0. However, it is simple 
to see that, up to a up to a field redefinition under the shift symmetry — > + 2ir, 
any regular solution of the differential equation (3.6) must obey this condition auto- 
matically. Asymptotically, as r — > 0, behaves as 

a „ Cr A - + Dr A + + ^- 
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with 



The requirement that A± is real gives rise to the Breitenlohner-Freedman bound 
m R L 2 > —1/4 for the bosonized scalar. Presumably after bosonization, m 2 R is pos- 
itive. However, it is tempting to speculate that bosonizing the fermions with the 
standard boundary condition (2.2) with — 1/2 < m < maps to a negative m 2 R above 
the BF-bound. If — | < m R L 2 < 0, then both fall-offs of are consistent with the 
= boundary condition, but we can specify a unique solution by demanding that 
the leading fall-off of is zero. 

Fractional Filling in the Hard Wall 

In Section (3.1), we imposed the boundary condition = 27m, with n G Z at the hard 
wall. However, as shown in [24], there is actually a one-parameter family of boundary 
conditions allowed at the hard wall. These arise by acting with a chiral symmetry, 
ipi = —e l9 ipR at the boundary. The effect of this chiral transformation on the boson 
is simply to shift the boundary condition to = 2irn + 9. These boundary conditions 
therefore impose that only fractionally filled Landau levels are allowed. 



C. Appendix: Field configurations in the near- horizon region 

In Section 3.2, we argued that regularity requires that the scalar field must attain an 
extremum of the effective potential Vm(4>) at the black hole horizon. 

2 

V IR {<f>) = m 2 R L 2 {\ - cos0) + %0 - O ) 2 

'"s 

When m 2 R L 2 > r^/r 2 , the effective potential has multiple extrema. To understand this 
situation further, we shall perform a detailed analysis of the field configurations in the 
near horizon AdS2 x M 2 region. Converting to the AdS2 radial coordinate 



c 



4 



Q(r h - r) 

we obtain the following equation of motion, valid asymptotically near the horizon. 

^ 9 , m 2 R L 2 r 2 
d 2 A = — ^— sin< 



g C 2 ^^^2 
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Note that this equation of motion respects the A0IS2 x IR 2 scaling symmetry ( — > X(, 
t — > Xt, B — > B. (The extra factor of 6 comes from the discrepancy between the AdS4 
and AdS2 radii). Exploiting the scaling symmetry, we cut off the AdS2 region at £ = 1 
and look for solutions in ( € [1, 00), satisfying the boundary condition <p(( = 1) = 0; 
other solutions are equivalent to these by scaling. 

The energy of the solution is captured by the effective Hamiltonian, 

/CO / 2 T 2 2 

dC + ^ (1 - cos0) + {(f) ~ 00)2 

All solutions to this equation of motion tend to a minimum of the effective potential 
but not necessarily to the global minimum. The reason is because the red-shift factor 
l/( 2 means that one can only gain a finite amount of energy by dropping into a lower 
minimum of the potential and it is possible for the gradient cost to outweigh the 
potential gain in making such a change. To illustrate this, suppose that <p\ and 02 are 
consecutive minima of Vtr(0), with 2 ~ 0i + 2n and V/#(0i) > V/r(0 2 ). Although 
ViR is lower at 02, it is not necessarily true that a solution that tends to 02 at the 
horizon is lower in energy than a solution tending to 0i; the field will save electrostatic 
energy by tending to 2 , but it costs energy to create the extra kink from X to 2 . 
Using standard flat space methods to estimate the energy of the kink, one finds that it 
is energetically favourable to create the extra kink to reach 2 if 

^ < ^2 - 0o) 2 - (0 2 - 0o) 2 ] (C.l) 

This is related to the pair production bound for fermions in AdS2 [46]. If one ignores 
the backreaction of the fermions, then it is only viable to place a kink at position r 
if the electrostatic gain ~ E/r is larger than than the effective mass of the fermion 
~ mft/r. Therefore pair production is only viable in AdS2 if one has E > m^. 

A typical situation is shown in Figure 8. When the filling fraction is in the vicinity of 
0o/27r ~ 4.83, Vj\r(0) has two minima at 0i/27r ~ 0.7, 02/27T ~ —0.1, with Vir{§\) > 
^tr(02)- For 0o/27r < 4.833, the minimum at 0i is preferred because the electrostatic 
energy saved by reaching 02 is not enough to compensate for the energy cost of creating 
the extra kink. Notice also that even though 0(£ — > 00) changes discontinuously as O 
crosses the critical value 4.833, the solutions 0(C; 0o) vary continuously for ( £ [l,oo) 
as 0o cross this threshold. As 0o increases, the position of the kink moves closer and 
closer to the horizon, and the kink disappears behind the horizon when O exceeds the 
threshold. 
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Figure 8: The electric field, E = r 2 F r t/ir = (</>o — 0)/2vr plotted against logarithm of the 
AdS 2 coordinate, In (, for filling fractions (Po/2tt = 4.9, 4.84, 4.834, 4.833, 4.82. The parameters 
are rn 2 R L 2 = 30 and r 2 Jr 2 = 6. 
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